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1. Abstract 2. GUI trajectory & Gravitational Waves 


We study binary black hole (BH) mergers in the extreme mass-ratio The merger scenario can be qualitatively divided into three stages: 
limit. ‘he energy, angular momentum and linear momentum of 
the remnant BH are determined based on the notion that toward the 
merger, small mass-ratio binary systems follow a geodesic universal in- 
fall (GUI) trajectory. Therefore, our analysis is performed directly in (ii) Universal plunge, where the infall path tends to the GUI trajectory. 
the test-particle limit by solving the Regge-Wheeler-Zerilli (RWZ) wave (iii) Quasinormal modes ringdown. 


equation with a source that moves along a geodesic. ‘This formalism cap- . . 

tures well the final inspiral stages of small mass-ratio binaries, and thus First wre solve aialy tically (he equations OF MOMOr (apy ute Ut © the 
’ Pp VAS aes Ce ie ae eodesic infall of a test particle, initially in a circular orbit at the ISCO 

provides a straightforward universal description In a region inaccessible to R= 6M). The GUL is universal in the sense that any small mass-ratio 

numerical relativity simulations. binary, that evolves due to GWs emission, tends to this infall trajectory 

after it crosses the ISCO. ‘Then, we calculate the GWs emitted by a test 

particle that moves along this geodesic by numerically solving the RWZ 


32. A Numerical Note - Extrapolation to Zt equation|1, 2}. The latter is a wave equation with a potential, induced by 
oa the spacetime curvature, and a source term, derived from the stress-energy 


The GWs flux is formally evalu- tensor: 
2,/,(A) 2 4/,(A) Ce Cy aera) 
ated at null infinity, Z*. Yet, we Ot Pe m — Ore Perm + Ve Vem = Sejm 


— —-—-R=250M 


(i) Quasicircular inspiral, during which the orbit evolves by the emission 
of gravitational waves (GWSs). 


can numerically extract them only 
at finite distances. We overcome 
this obstacle by analytically calcu- 
lating the GWs propagation from a 
given distance to Z* by solving the 
homogeneous RWZ equation in the 
frequency domain. 

As can be seen in Figure 2, after 
the extrapolation to Z*, the align- | 


ment of the (,W amplitude that Figure 2: GWs extracted at R = 250M wqo00 4200 4400 4600. 4800 000 6200 6400 


was numerically extracted at dif- (blue dashed line), 750M (red dash-dotted : ana 
ae line), 1500M (yellow dotted line), and their 
ferent radii improves by more than extrapolation to Z+ (solid lines, respective Figure 1: The left panel presents the GUI trajectory of a test particle that inspirals 


an order of magnitude. colors). from the ISCO into the horizon (black dashed line). The right panel presents the 
corresponding GWs emission. 


4. Results: Recoil Velocity 


Using the above GWs calculation, we evaluate the linear momentum flux and the corresponding recoil velocity The test particle GUI waveform, 
of the remnant BH. We calculate numerically the contribution up to € = 10 multipole. Assuming that the shown in the right panel of Figure 
contributions of higher multipoles decay exponentially, as can be seen in Figure 3, we estimate the total recoil 1, has a linear dependence on the 


5. Waveform 


velocity, which has a quadratic dependence on the symmetric mass ratio, v, and is given by: mass ratio and hence it can be 
easily scaled to capture the late 

V/c & 0.0467*. stages of realistic small mass-ratio 

This result is greater than the known value in the literature[3, 4] by approximately 4%, mostly as a result of the binary mergers. ‘The dephasing 
high multipole contributions. time of the GUI and the small 


mass-ratio signal scales as: 


tap e< py— 3/10. 
0.046 Thus, for example, given a merger 
of a stellar mass object and a 
0.044 super-massive BH, with a mass of 
ees 10° — 10°Mo, the GUI waveform 
: faithfully describes between tens to 
0.04 hundreds of the final GWs cycles. 
N 
so 6. Final Mass & Spin 
= 0.036 The final mass and spin of the 
0.034 remnant BH are given by: 
0.032 
e Num. result ). 
0.03 ae 
_ VV © 0.0467 The first order term is derived 
0.028 analytically from the energy and 
. : : : . : : 2 uw angular momentum of the test 
£ particle at the ISCO: 
Figure 3: The final recoil velocity, up to a given multipole @. The total velocity converges to V/v* = 0.0467c (black dash-dotted Ersco/- = Je meio ey (a Ie 


line) and corresponds to an exponential decay trend (red dashed line) The next-order term corresponds 


to the GWs emission after the 
ISCO crossing. Its scaling can be 


eferences deduced based on the scaling of 


* and 


the emitted fluxes, F «x v 
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